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Measurement of the gravitational constant G in space (Project SEE): 
sensitivity to orbital parameters and space charge effect 

A.D. Alexeev^, K.A. Bronnikov^, N.I. Kolosnitsyn^, M.Yu. Konstantinov'^, 
V.N. Melnikov^t and A.J. Sanders^ 

f Russian Gravitational Society, 3-1 M. Ulyanovoy St., Moscow 117313, Russia 

t Dept. of Physics and Astronomy, University of Tennessee, Knoxville, TN 37996-1200, USA 

We describe some new estimates concerning the recently proposed SEE (Satellite Energy Exchange) experiment for 
measuring the gravitational interaction parameters in space. The experiment entails precision tracking of the relative 
motion of two test bodies (a heavy "Shepherd", and a light "Particle") on board a drag-free space capsule. The 
new estimates include (i) the sensitivity of Particle trajectories and G measurement to the Shepherd quadrupole 
moment uncertainties; (ii) the measurement errors of G and the strength of a putative Yukawa-type force whose 
range parameter A may be either of the order of a few meters or close to the Earth radius; (iii) a possible effect of 
the Van Allen radiation belts on the SEE experiment due to test body electric charging. The main conclusions are 
that (i) the SEE concept may allow one to measure G with an uncertainty smaller than 10~^ and a progress up to 
2 orders of magnitude is possible in the assessment of the hypothetic Yukawa forces and (ii) van Allen charging of 
test bodies is a problem of importance but it may be solved by the existing methods. 



1. Introduction 

The SEE (Satellite Energy Exchange) concept of a 
space-based gravitational experiment was suggested in 
the early 90s ^ and was aimed at precisely measuring 
the gravitational interaction parameters: the gravita- 
tional constant G , possible violations of the equivalence 
principle measured by the Eotvos parameter 77, time 
variations of G , and hypothetical non-Newtonian grav- 
itational forces (parametrized by the Yukawa strength 
a and range A). Such tests are intended to overcome 
the limitations of the current methods of ground-based 
experimentation and observation of astronomical phe- 
nomena and to fill gaps left by them. The significance 
of new measurements is quite evident since nearly all 
modified theories of gravity and unified theories pre- 
dict some violations of the Equivalence Principle (EP), 
either by deviations from the Newtonian law (inverse- 
square-law, ISL) or by composition-dependent (CD) 
gravity accelerations, due to the appearance of new 
possible massive particles (partners); time variations of 
G (G-dot) are also generally predicted ^, ||]. 

Since gravitational forces are so very small, precision- 
measurement techniques have been at the core of ter- 
restrial gravity research for two centuries. However, 
evidence is increasingly accumulating which indicates 
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that terrestrial methods have plateaued in accuracy and 
are unlikely to achieve significant accuracy gains in the 
future. Ref. [Q describes the contemporary situation 
relative to the measurements of G and G-dot. 

There is considerable evidence that the uncertainty 
in G has plateaued at about 100 ppm |[ 0, Most 
of the stated errors in the recent experiments are of the 
order 100 ppm. Moreover, the scatter (1-sigma) about 
the mean is about 140 ppm. This conclusion does not 
include the recent paper by Gundlach and Merkowitz 
|9j who report an error of about 14 ppm in G measure- 
ment using a dynamically driven torsion balance. The 
reported value of G agrees with the previous measure- 
ments within their uncertainties. This work is unique 
for the moment and needs confirmation by other exper- 
imental groups using alternative methods. 

It might seem that the problems of terrestrial appa- 
ratus must inexorably yield to new technologies — that 
the promise of ever increasing sensitivities would also 
lead to ever improving accuracy. However, this may not 
be true, since it is the intrinsically weak nature of the 
force and the resulting systematic errors which arise in 
its isolation and measurement that limit the ultimate 
attainable accuracy in terrestrial experiments Q . 

The EP may be tested by searching for either vio- 
lations of the inverse-square law (ISL) or composition- 
dependent (CD) effects in gravitational free fall. 

In the watershed year of 1986, Fischbach et al. star- 
tled the physics community by showing that Eotvos's 
famous turn-of-the-century experiment is much less de- 
cisive as a null result than was generally believed | [To| , 
p| . Prior to this time, experiments by Dicke |l^ and 
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Braginsky [[l3[ had demonstrated the universahty of 
free faU (UFF) to very high accuracy with respect to 
several metals falling in the gravitational field of the 
Sun (the Eotvos parameter 77 was ultimately found to 
be smaller than 10^^^). The interpretation of these 
results at the time was that they validated UFF. 

Since 1986 it has become customary to parametrize 
possible apparent EP violations as if due to a Yukawa 
particle with a Compton wavelength A . This approach 
unites both ISL and CD effects very naturally, while 
the parameter values in the Yukawa potential suggest 
which experimental conditions are required to detect 
the new interaction. 

Following the conjecture of Fischbach et al., ISL and 
CD tests were undertaken by many investigators. Al- 
though a number of anomalies were initially reported, 
nearly all of these were eventually explained in terms 
of overlooked systematic errors or extreme sensitivity 
to models, while most investigators obtained null re- 
sults. However, a positive result for a deviation from 
the Newtonian law (ISL) was obtained (and interpreted 
in terms of a Yukawa-type potential) in the range of 20 
to 500 m by Achilli and colleagues iQ; this needs to 
be verified in other independent experiments. 

For reviews of terrestrial searches for non- Newtonian 
gravity, see ^ 0. 

The idea of the SEE method is to study the relative 
motion of two bodies on board a drag-free Earth satel- 
lite using horseshoe-type trajectories, previously well- 
known in planetary satellite astronomy: if the lighter 
body (the "Particle" ) is moving along a lower orbit than 
the heavier one (the "Shepherd" ) and approaching from 
behind, then the Particle almost overtakes the Shep- 
herd, but it gains energy due to their gravitational in- 
teraction, passes therefore to a higher orbit and begins 
to lag behind. The interaction phase can be studied 
within a drag-free capsule (a cylinder up to 20 m long, 
about 1 m in diameter) where the Particle can loiter 
as long as 10^ seconds. It was claimed that the SEE 
method exceeded in accuracy all other suggestions, at 
least with respect to G and a for A of the order of 
meters. Some design features were considered, making 
it possible to reduce various sources of error to a neg- 
ligible level. It was concluded, in particular, that the 
most favorable orbits are the sun-synchronous, contin- 
uous sunlight orbits situated at altitudes between 1390 
and 3330 km 0. 

Since the origination of the SEE concept, the de- 
velopment has focused on critical analyses of orbital 
parameters and satellite performance and the assess- 
ment of critical hardware requirements. All indications 
from this work are that the SEE concept is feasible and 
practicable [ p8[ . 

At the present stage one can assert that, although 
space is a challenging environment for research, the in- 
herent quietness of space can be exploited to make very 



accurate determinations of G and other gravitational 
parameters, providing that care is taken to understand 
the many physical phenomena in space which have the 
potential to vitiate accuracy. A distinctive feature of a 
SEE mission is its capability to perform such determi- 
nations simultaneously on multiple parameters, making 
it one of the most promising proposals. 

To be more specific, let us enumerate the suggested 
SEE tests and measurements and show their expected 
accuracy as currently estimated: 



Test/measurement 

EP/ISL at a few metres 
EP/CD at a few metres 
EP/ISL at - i?e 
G 

G/G 



Expected accuracy 
2-10-^ 

< lO-'^ (a < 10-4) 

< 10-1" 
3.3-10-^ 



< 10 



-13 



m one year 



The last estimate is only tentative; the subject is 
under study. 

This paper presents a description of the main ideas 
of the SEE experiment and some new evaluations con- 
cerning the opportunities of the SEE concept and its 
yet-unresolved difficulties. In Sec. 2, on the basis of 
computer simulations of Particle trajectories, we de- 
scribe some general properties of Particle trajectories in 
short (Particle-Shepherd separation from 2 to 5 meters) 
capsules, their sensitivity to the value of the parame- 
ters of the gravitational interaction, and we estimate 
the requirements to the Shepherd quadrupole moment 
uncertainty. Sec. 3 shows the results of simulations of 
the measurement procedure itself, which enables us to 
estimate the possible measurement accuracy with re- 
spect to G and a for A of the order of either meters 
or the Earth's radius. Sec. 4 discusses a spurious effect 
of test body electric charging when the satellite orbit 
passes through the Van Allen radiation belts, rich in 
high-energy protons. Sec. 5 is a conclusion. 

In what follows, the term "orbit" applies to satel- 
lite (or Shepherd) motion around the Earth, while the 
words "trajectory" or "path" apply to Particle motion 
with respect to the Shepherd inside the capsule. 

2. Simulations of Particle trajectories 

In the previous studies of the SEE project it was as- 
sumed that the capsule was about 20 m long and the 
initial Shepherd-Particle separation xq along the cap- 
sule axis was as great as 18 m; some estimations were 
also made for 5 m < .xo < 10 m. The Shepherd mass 
was taken to be M = 500 kg and the Particle mass 
rn ~ 0.1 kg. The present study retains these values, 
except the cases noted below. 
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In this section we will describe some characteristic features of Particle trajectories in a short capsules (Particle- 
Shepherd separation 2 m < a;o < 5 m) with the Shepherd mass reduced to M = 200 kg. Our goal is to determine 
the properties of trajectories in the case considered and to determine the sensitivity of trajectories to the uncertainty 
of orbit radius, the value of the Newtonian gravitational constant and of the Shepherd quadrupole moment J2 • As 
in our previous studies, the capsule diameter is presumed to be 1 m. 

The reason for considering the short capsule, the Shepherd with reduced mass and the quadrupole moment 
uncertainty is economical and technological in origin. Namely, it is hard to produce a spherically symmetric 
Shepherd to the required accuracy. To avoid the inclusion of SJ2 in the set of parameters to be determined in the 
experiment, it is useful to know which values of SJ2 will be negligible, since the growth of the number of parameters 
leads to serious problems in data processing. 



2.1. Equations of motion and initial data 



For simplicity, we assume that the relative motion of the test bodies inside the capsule occurs in the satellite orbital 
plane. This simplification is purely technical, since, as was found in our previous studies, the three-dimensional 
nature of the Particle motion does not change the main estimates. 

The reduced Lagrangian of the Particle motion in the considered case is 



G 



M^m Mm 



G^{l + J.(^)V.(cos^)}(l 



(1) 



where (R, ip) are the Earth-centered polar coordinates of the Shepherd in the orbital plane; r = ^ {R-VyY' -Vx^ and 
V' are the Earth-centered polar coordinates of the Particle; x and y are the Shepherd-centered Particle coordinates, 
where x is the "horizontal" one, i.e., along the orbit and simultaneously along the capsule and y is the "vertical" 
one, along the Earth-Shepherd radius vector; s = sjx^ ^y^ is the Particle-Shepherd separation; , M and m 
are the Earth, Shepherd and Particle masses, respectively; J2 is the quadrupole moment of the Shepherd, Tg is its 
radius and Pi is the Legcndre polynomial 



P2(cos6l) = 



3 cos^ d - 1 



where 6 is the angle between the line connecting the centres of test bodies and the Shepherd equatorial plane. It 

is easy to see that if the Shepherd symmetry axis is in its orbital plane, then = 9q = — arctan(j//.x) + cp. If the 
symmetry axis of the Shepherd is orthogonal to its orbital plane, then 6* = . In general, if x is the angle between 
the Shepherd symmetry axis and its orbital plane, then = 60 cos x ■ Hence the influence of J2 on the Particle 
motion is minimum if the Shepherd symmetry axis lies in its orbital plane and is maximum if they are mutually 
orthogonal. 

For simplicity (and taking into account the corresponding estimate) we neglect the influence of the Particle on 

the Shepherd, so the Shepherd trajectory is considered as given. Then, varying the above Lagrangian with respect 
to X and y , taking into account that M ^ m and i? ^ s , we arrive at the following equations of Particle motion 
with respect to the Shepherd: 



fx 



2yip -\-x{ip 



2R(py CM 



a;|l + J2 (y)'p2(cos^)| 



CM 



— ax- 



+ ^^J° J2 (l + ae-^/^^ X [a;(l-|-3cos26i) -F3t/sin26lcosx] ; 



(2) 



^ = -2xip + {R + y)-{ip' - 



+ 



2Ripx GM 



-y 



1 + J2 (^)V2(COS0) 



ay- 



GM (\ 



^ ,i){l + J2(^)V2(cOS.)^ 

^^^J2 (l-hae-*/-^) X [3a: sin 26* cos x + 2/(1 -3 COS 61)] 



(3) 



where M = M + m. 
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Two kinds of initial conditions for Eqs. (g) and (j^) 
were used during the simulations. First, we used the so- 
called "standard" initial conditions, taking the Particle 
velocity components i(0) and y{0) corresponding to 
its unperturbed (i.e., without the AI — rn interaction) 
orbital motion distinguished from the Shepherd's orbit 
only by its radius (for circular orbits) or semimajor axis 
(for elliptic orbits). Assuming that the Particle motion 
begins right at the moment when the Shepherd passes 
its perigee, these conditions have the form 



x(0) 
i(0) 



2^0, y{0) 

2(1 -e)2' 



yo, 

m 



ujexQ 
'e'(l -e) 



(4) 



where — GM^/Rq, Rq is the Shepherd orbital ra- 
dius (at the perigee), e is the orbital eccentricity and 

For clearness, the relations (^) are written in the 
linear approximation in the variables x and y . Higher- 
order approximations were used in the simulation pro- 
cess as well. 

The second kind of initial conditions corresponds to 
small variations of initial velocities with respect to their 
"standard" values. 

The set of equations (||)-(||) was solved numerically 
using the software developed previously ||l8| to analyze 
the SEE project. 

On the basis of numerical solutions of Eqs. (||) and 
(^, we considered two types of Particle trajectories, 
corresponding to different choices of the initial data: 
(i) approximately U-shaped ones and (ii) cycloidal ones, 
containing loops (see more details on the trajectories in 
[§ |I|]), for orbital altitudes iJorb = 500 , 1500 and 3000 
km. 

2.2. General characteristics of trajectories in 
a short capsule 

Reduction of the Shepherd mass leads to the reduction 
of distance between libration points Li and L2 . As a 
result, the region where the horseshoe orbits exist is re- 
duced as well, and the turning points for the horseshoe 
orbits become closer to the Shepherd. Thus, for a Shep- 
herd mass of M — 200 kg, the positions of the turning 
points of the horseshoe orbits starting, for instance, at 
xo = 18 m and 0.1 m< |yo| l£ 0.3 m are placed in the 
region a; < 5 m for all orbits whose altitude H is in 
interval 500 km < H < 3000 km. The usage of short 
trajectories, which start in the region xq < 5 m, and 
especially extremely short trajectories with xq = 2 m, 
leads to additional limitations on the initial conditions. 

Numerical simulation shows that to avoid Particle 
collisions with the Shepherd and the walls of the cap- 
sule for "smooth" trajectories, which correspond to the 
standard initial conditions, trajectories with \yo\ < 0.2 
m may be used for orbit altitudes Ho^h ^ 1500 km 



while for orbit altitude iJorb = 3000 km and xq > 3 m 
trajectories with |yo| = 25 cm may be used also. 

The usage of cycloidal trajectories, which appear 
when the absolute values of the Particle's initial velocity 
exceed its standard value, make it possible to avoid this 
limitation. 

2.3. Trajectory sensitivity with respect to the 
Newtonian gravitational constant 

The influence of the Newtonian gravitational constant 
G on the Particle motion in the case considered may 
be investigated by considering how a small perturbation 
SG from the "standard" value of G changes the Particle 
trajectories. Such a perturbation is characterized by 
the displacement of a perturbed trajectory with respect 
to an unperturbed one, that is, 

Sr^ fs{t) - ro{t) 

where fs and fg denote the Particle radius vector for 
perturbed and unperturbed motion, respectively. For 
rather long trajectories, instead of the full displacement 
6f, the displacement along the x axis ( Sx ) may be con- 
sidered, because a numerical simulation shows that a 
displacement along the y axis is one order of magni- 
tude smaller than a displacement along the 2; -axis. 

U-shaped Particle trajectories were considered for 
circular Shepherd orbits in the following range of pa- 
rameters and initial conditions: orbit altitudes i?orb = 
500, 1500 and 3000 km; the initial position changes 
in the range 2 m < xq < 5 m, —20 cm < yo ^ —5 
cm. For comparison with the case where M — 500 kg, 
trajectories with = 18 m were considered also. It 
was found that the perturbation SG/G — 10^^ leads to 
displacement of smooth trajectories in the range from 
0.634-10~® mto 3.34-10"^ m with a minimum displace- 
ment achieved for yo — 0.05 m, = 2 m and i?orb = 
3000 km while a maximal displacement is achieved for 
yo — 0.05 m, = 5 m and i/orb = 500 km. 

The use of cycloidal trajectories increases the influ- 
ence of 6G on the Particle motion. 

2.4. Trajectory sensitivity to the orbit radius 
uncertainty 

This sensitivity is characterized by the displacement 
6x , induced by a small perturbation (or uncertainty) 
6h of the orbital altitude |^ -ffoib • 

It was found that for Sh ~ 1 cm, the orbital altitude 
^orb = 1500 km and \yo\ = 20 cm, the maximum 
value of 6x increases from 3.222-10^^ mat xo = 2 
m to 9.387 • 10"*^ m at = 5 m. For H^rh = 500 
km these values must be multiplied by a factor of 2 



^As above, we consider the displacement 5x instead of Sr 
because 5x provides the main contribution to Sr , while Sy is 
much smaller. 



Measurement of the gravitational constant in space (Project SEE) 



5 



Table 1. Estimates of 5G/G in ppm for 6J2 = 10"^ , 
when the symmetry axis of the Shepherd is orthogonal 
(x = ^f^) to its orbital plane. The second line shows 



yo, 

cm 




2 m 


3 m 


4 m 


5 m 


-5 


1.568 


5.13 


2.41 


0.138 


-10 


3.181 


1.41 


0.876 


0.649 


-15 


0.857 


6.08 


5.057 


4.486 


-20 


60.28 


44.31 


31.91 


28.04 



(approximately) and for iforb = 3000 km by a factor 
of 1/2. 

It was also found that the dependence Sx{5h) is, to 
a good accuracy, linear: in particular, for 6h = 1 m, 
the orbital altitude iJorb = 1500 km, a;o = 5 m and 
\yo\ = 20 cm, 6x = 9.387 • lO"*' m as expected. 

The use of cycloidal trajectories reduces the depen- 
dence of trajectories on the orbit altitude: for instance, 
for iJorb = 1500 km, xo = 5 m and \yo\ = 20 cm, the 
uncertainty Sh ^ 1 cm leads to Sx ~ 1.455 • 10^® m. 

It is necessary to point out that the influence of 
the orbital altitude on the Particle trajectory for the 
reduced Shepherd mass (M = 200 kg) is greater than 
in the case M = 500 kg. For instance, in the case 
M = 500 kg, Hoib = 1500 km, Sh = 1 m, xq = 5 m 
and |t/o| =20 cm, the trajectory displacement becomes 
6x = 4.6593 • 10-8 m. 

These estimates lead to certain restrictions on the 
reasonable precision of Particle trajectory measure- 
ments. 

2.5. Trajectory sensitivity to the Shepherd's 
quadrupole moment uncertainty 

A maximum effect of the Shepherd quadrupole moment 
J2 on the Particle motion is realized in the case when 
the Shepherd axis is orthogonal to its orbital plane. 
The influence of SJ2 on the accuracy of G measure- 
ment may be estimated as follows. Let some value 
of 6J2 produce the trajectory displacement \Sr\ < Slj 
while the variation SGq of G with the same initial con- 
ditions gives the trajectory displacement \6r\ < SIq ■ 
Then, keeping in mind the linear dependence of trajec- 
tory displacements on SJ2 and 6G , the accuracy of G 
measurement under the uncertainty 5J2 may be esti- 
mated as 



SG 
~G 



< 



1 -I 

SIg 



6G0 
G ' 



Using this inequality and the results of trajec- 
tory simulations, we obtain the following estimates for 
U-shaped Particle trajectories in circular orbits with 
Horh = 1500 km: 



The uncertainties 5J2 ^ 10^'"' do not create a sub- 
stantial error in G for most of the trajectories. Only 
trajectories with |yo| = 0.2 m require 5J2 ^ 10^^ be- 
cause of the growth of the sinusoidal component for 
these trajectories. 

3. Simulations of experimental 
procedures 

This section describes the results of numerical simula- 
tions of the whole measurement procedures aimed at 
obtaining the sought-after gravitational interaction pa- 
rameters. These simulations assumed a Shepherd mass 
of M = 500 kg, a circular orbit with i?orb = 1500 km 
with a spherical gravitational potential for the Earth, 
and a Particle mass of 100 g. Where relevant, it is 
assumed that both the Shepherd and the Particle are 
made of tungsten. Identical compositions for them are 
assumed for simplicity since this work is performed for 
estimation purposes only. 

3.1. Simulations of an experiment for 
measuring G 

The constant G is determined from the best fitting 
condition between the "theoretical" Particle trajecto- 
ries (f^^iU) = f/'^), calculated by Eqs. (2) and (3) for 
J2 = and the "empirical" (ri) Particle trajectories 
near the Shepherd. The fitting quality is evaluated by 
minimizing a functional that characterizes a "distance" 
between the trajectories. We have considered the fol- 
lowing functional for such "distances" : 



N r 



{x,-xf^r + {yi-yfr 



(5) 



i=l 

N 

5. = ^(x,-^f)2, Sy = J2iyi-yfr, (6) 



i=l 



The theoretical trajectory depends on the gravi- 
tational constant G, on the initial coordinates Xo,yo 
and on the initial velocities Vxo,Vyo- To estimate G, 
one chooses the value for which a "distance" functional 
in the space of the five variables {G, xo,yo,'Vxo,Vyo) 
reaches its minimum. 

We carried out a numerical simulation of the SEE 
experiment and estimated 6G for a given coordinate 
measurement error (cr = I-IO"^ m). As "empirical" 
trajectories, we took computed trajectories, with spec- 
ified values of the above five variables, where Gaussian 
noise was introduced from a random number genera- 
tor. Independent "empirical" trajectories were created 
by non-intersecting random number sequences. The 
functional was minimized using the gradient descent 
method and the consecutive descent method. The start- 
ing value of the "vertical" (along the Earth's radius) 
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Figure 1: Uncertainties SG estimated by the gradient de- 
scent (-Rgrad) and consecutive descent (-Rs) methods. 



coordinate, yo, was taken to be 0.25 m, while the hori- 
zontal one, xq , varied between 2 and 18 m. Fig 1 shows 
the dependence of the errors SG/G — Rgrad, obtained 
by the gradient descent method and SG/G — Rg, ob- 
tained by the consecutive descent method. All the er- 
rors are estimated by confidence intervals corresponding 
to a confidence of 0.95. The mean values of these errors 
are as follows: 



R, 



grad 



4.69•10-^ 



Rs 



5.24-10-8. 



Thus the errors estimated by the gradient and consec- 
utive descent methods are close to each other and are 
about an order of magnitude smaller than the error 
from one-trajectory data. It was found that the sim- 
ulation results depend strongly on the random num- 
ber generator, so that ordinary generators are not per- 
fect: the generated random number sets do not obey 
the Gaussian law. 

The use of truncated functionals like (2) has shown 
that a functional incorporating the more informative 
"horizontal" coordinate x leads to estimates close to 
those obtained from the total functional, whereas the 
use of y alone substantially decreases the sensitivity. 
Therefore in practice, to determine G, it is sufficient to 
measure only one of the two coordinates, viz. x. 

Since the "empirical" trajectory is built on the basis 
of a computed one, with a known value of the gravita- 
tional constant Go , it appears possible to estimate a 
possible systematic error inherent in the data process- 
ing method. The latter has turned out to be in most 
cases much smaller than the random error. This result 
shows the correctness of the methods used. 

As is evident from the results, the best accuracy is 
achieved at values of xq ( ~ the capsule size) about 4-5 
meters. 
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Figure 2: Random "velocity" and "coordinate" uncertain- 
ties. 



3.2. Simulation of measurement procedures 
for estimation of G using velocity data 

When one can precisely measure the Particle velocity, 
the gravitational constant can be estimated from the 
velocity data alone. A computer program was created 
for simulating an experiment determining the gravi- 
tational constant G by velocity data. These simula- 
tions assumed a Shepherd mass of 500 kg, a circular 
orbit with iJorb = 1500 km (spherical Earth's poten- 
tial), a Particle mass of 100 g, and a velocity error 
8vx — 1 • lO^'* m/s. The constant G is estimated from 
the best fitting condition between the theoretical Parti- 
cle velocity {y^'^iti) = vf^) and the empirical one. The 
fitting quality is evaluated by minimizing a functional 
characterizing a distance between velocity trajectories. 
We used the following functional: 



V 



N 



The theoretical velocity depends on the gravita- 
tional constant G , on the initial coordinates xq , yo 
and the initial velocities w^oi^^yo- To estimate G, one 
chooses the value for which the velocity distance func- 
tional V reaches its minimum in the space of the five 
variable ( G, xq , yo , i^xo, ) 

As empirical trajectories, we took computed trajec- 
tories with specified values of the above five variables, 
where Gaussian noise was introduced from a random 
number generator. Independent empirical trajectories 
were created by non-intersecting random number se- 
quences. The functional V was minimized using the 
gradient descent method. The starting value of the yo 
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coordinate was taken to be 0.25 m, while the horizontal 
one, xo , varied between 2 and 6 m. The value of G was 
estimated from 11 trajectories at a significance level of 
0.95. The appropriate errors were estimated at confi- 
dence intervals and then related to G as shown by the 
upper curve in Fig. 2. The relative error obtained using 
the coordinate functional for a coordinate measurement 
error of cr = 1 • 10^^ m is also shown in the same figure. 
Computer-based simulation also allows one to estimate 
systematic errors. We see that the accuracy of estima- 
tion of G for a given set of measurement errors is the 
best for the coordinate functional. 

3.3. Equations of motion with Yukawa terms 

We will present the Particle equations of motion in 
the relevant approximation, including the contributions 
from hypothetical Yukawa forces, taking into account 
the finite size of the Yukawa field sources. 

Let the interaction potential for two elementary 
masses mi and m2 be described by the potential 



dV 



Yu 



G dmidm2 



-a e 



-r/\ 



(7) 



where r is the masses' separation, a and A are the 
strength parameter and the range of the Yukawa forces. 
Then for two massive bodies with the radii Ri and i?2 
after integration over their volumes we obtain pC| 



^Yu ^ Grnim2PiP2 ^^-r/x 



where 



Ri , Ri . , Ri 

—— cosh — smh — 

A A A 



(8) 



(9) 



When Ri/X <^ 1 , we have /3i « 1 . This may be the case 
when we consider the interaction between the Shepherd 
and the Particle at a distance of the order of a few 
meters. The radii of the Shepherd and the Particle are 
small: i?i « 18 cm for the Shepherd and i?2 ~ 1-1 
cm for the Particle. If the range A is of the order of 
the Earth radius, A ~ R^, we have = 1.10 and 
/3i.2 = 1 where the indices 1 and 2 label the Shepherd 
and the Particle, respectively. 

Equations of motion are obtained under the follow- 
ing assumptions. There are two Yukawa interactions 
with the parameters Aq and ao referring to the Earth- 
Shepherd and Earth-Particle interactions which are the 
same (due to the assumed identical composition for the 
Shepherd and the Particle), while A and a determine 
the Shepherd-Particle interaction. The equations of 
motion in the frame of reference of the Shepherd, with 
the same notations for x, y and s as used previously, 
are 

X + 2oJ^y + G{mi + ^2)^ — Sui^ — 



G(mi + m2)4a 1 + T ) ^"'^^ = 



y - 2ujx - 3uj^y + G(mi + 1712)-^ + ( y^ 



X 

"2" 



+ G{mi + m2)^a(^l + e^'/^ 

-tj2/3oaoe~''"i/^°?/ = (10) 



where u is the orbital frequency: 



GMff 



(11) 



We have neglected the terms quadratic in s/vqi times 
a or ao due to their manifestly small contributions. 



If we set ao = in Eqs. (10), we obtain the equa- 



tions used to describe only the Shepherd-Particle Yukawa 
interaction. Notice that the Yukawa terms are roughly 
proportional to the gradients of the corresponding New- 
tonian accelerations, namely, Grrii/s^ for the Shepherd- 



Particle interaction and GMg)/r\ 



for (say) the 



Earth-Shepherd interaction. In our case these quanti- 
ties are estimated as 



Gmi 



2.7-10-1" s-2 



and 



for s = 5 



8.16-10-^ s-^. 



(12) 



Thus, given the same strength parameter, the Earth's 
Yukawa force is three orders of magnitude greater than 
that between the Shepherd and the Particle. Therefore, 
one might expect some significant progress in an ISL 
test for A of the order of the Earth's radius. 

Eqs. ([lo| ) were used to simulate the measurement 
procedures. 

3.4. Sensitivity to Yukawa forces with A ^ 1 m 

In an experiment for finding the Yukawa interaction 
between the Shepherd and the Particle using the po- 
tential (^ with /?i_2 = 1, one computes two theoret- 
ical trajectories: the first ignoring the Yukawa forces 
(a;"(ii), y'^{ti)) and the second taking them into ac- 
count (^x^iti), y^iti)^ . These two computed curves 
are compared with the empirical trajectory using the 
functional Sk (fc = 0, a) according to (||) which may be 
considered as a dispersion characterizing a scatter of the 
"empirical" coordinates with respect to the fitting tra- 
jectory. This is true when the theoretical model is ade- 
quate to the real situation. In the case k ^ a the func- 
tional Sk = Sa has a distribution with 712 = 2N — 1 
degrees of freedom. With k = the parameter a is 
absent, therefore 5*0 is distributed according to the 
law with iVi = 2N degrees of freedom. Then their ra- 
tio So/Sa — Fn2,ni wiU be distributed according to the 
Fischer law [E^ with n2 and ni degrees of freedom. If 
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an experiment shows that, on a given significance level 
q, the relation (see Q) 



So/Sa ^ Fni,n2,q 



(13) 



is valid, one should conclude that the Yukawa force has 
been detected. An equality sign shows a minimum de- 
tectable force for the given significance level q . We have 
assumed q = 0.95 . The results of a sensitivity compu- 
tation for different values of the space parameter A are 
presented in Fig. 3 (curve 1). A maximum sensitivity 
of a — 2.1-10"^ has been observed for A = 1.25 m. 
This value is 3 to 4 orders of magnitude better than 
the sensitivity of terrestrial experiments in the same 
range. 

3.5. Sensitivity to Yukawa forces with A ^ 

To estimate the parameter ao in Eqs. (p^), computer 
simulations were carried out using the method as de- 
scribed above for a, based on the Fischer criterion for 
the significance level 0.95. The range parameter Aq 
varied from (l/32)i?0 to 32i?0 . Two trajectories with 
the initial Shepherd-Particle separations xq of 2 and 5 
m were calculated. In both cases the impact param- 
eter uq was chosen to be 0.25 m. We used Eqs. ( [lO| ) 
with a ~ 0, i.e., excluding the non-Newtonian inter- 
action between the Shepherd and the Particle. As is 



evident from Eqs. (10), the Particle trajectory depends 
on the ratio roi/Ao in the product (roi/ Xq) e~^°^^^° . 
This quantity reaches its maximum at Aq = roi/2 . Our 
calculations have confirmed that a maximum sensitiv- 
ity of the SEE method (3.4-10"'^ for xq = 5 m) is 
indeed observed at this value of Aq . This is about 
an order of magnitude better than the estimates ob- 
tained by other methods. Hopefully this estimate may 
be further improved by about an order of magnitude 
by optimization of the orbital parameters. However, 
there is a factor which can, to a certain extent, spoil 
these results, namely, the uncertainty in the parameter 
w which, in this calculation, was assumed to be known 
precisely. 

The simulation results are shown in Fig. 3 (curve 2) 
for a trajectory with an initial Shepherd-Particle sep- 
aration of 5 metres. See the next section for further 
discussion of this figure. 

3.6. Precession of the Shepherd orbit and a 
test of the Inverse Square Law at 
distances of the order of i?® 

We have shown previously that the SEE experiment al- 
lows one to test the inverse square law at distances on 
the order of 1 m (amin ~ 2 • 10^'') and at distances on 
the order of a half of an orbit radius (amin ~ 3.4- 10^^) 
[ p5[ . Another test can be done using spacecraft preces- 
sion data. As known, in two-body problems an orbit is 



closed for only two potentials. They are (1) the New- 
tonian potential, U ^ l/r , and (2) U ^ . In other 
cases the orbit is not closed and a pericenter precession 
is observed. In particular, any deviation from the New- 
tonian law entails a precession of an orbit due to the 
Yukawa interaction 



t/-=^™ae-/^ 



(14) 



the Shepherd orbit exhibits a precession. In the general 
case the precession magnitude due to a small perturba- 
tion, described by a potential SU , is given by (see [1]) 



(15) 



Integration is done over a non-perturbed trajectory. 
Here m is the Shepherd's mass, m' is the mass of 
a central body (the Earth), M = mr'^ip is an inte- 
gral of motion (the angular momentum), and 5U = 
a{Gmm') e~'^^^ . The non-perturbed trajectory is de- 
scribed by the expressions 

p 

P 



1 + e cos ip 

2EAP 



ni^Gmm') 



j-p = a{l~e ). 



(16) 



After a standard algebraic computation we obtain 



S(p = a— 



e Jo (1 + e cos ipy 



1 



X < - [2e + (1 + e^) cos ip] - (e + cos ip) ^ dip, (17) 



where 



— (18) 
A A 1 + e cos ip ' 

Using Eq. (8) and data with the error Sp for the SEE 
Satellite, we calculated the curves a{X) , which deter- 
mine the border on the a — A plane between two do- 
mains, where the Yukawa interaction (a new long-range 
force) is forbidden by experiment and where it is not. 
The sensitivities to Yukawa interactions are shown in 
Fig. 3 as the curve 4 (see Q ) for the parameter A in 
the range from 1 • 10^ m to 1 • 10^^ m. The curve 
3 was calculated for the SEE satellite with an eccen- 
tricity e = 0.01 and the precession error is equal to 
Sp, = 0.1"/y. 

One can conclude that the inverse square law may 
be tested with a sensitivity ofQ!^6.3T0^^^ forA~ 
3.9 TO^ m (half the orbital radius). 

4. A possible effect of the Earth's 
radiation belt 

Charged particles, penetrating into the SEE capsule 
from space and captured by the test bodies, create elec- 
trostatic forces that could substantially distort the ex- 
perimental results. Among the sources of such particles 
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Figure 3: The SEE method sensitivity to Yukawa forces with the range parameter Ao of the order meters (1) of the Earth's 
radius R(g using trajectory measurements (2) and orbit precession (3). The limitations on Yukawa's forces parameters from 

(4)- 



one should mention (i) cosmic-ray showers, (ii) solar 
flares and (iii) the Earth's radiation belts (Van Allen 
belts). The effect of cosmic-ray showers was estimated 
in Ref. and shown to be negligible. Solar flares are 
more or less rare events and, although they create very 
significant charged particle fluxes, sometimes even ex- 
ceeding those in the most dense regions of the radiation 
belts, one can assume that the SEE measurements (ex- 
cept those of G) are stopped for the period of an intense 
flare. On the contrary, the effect of the Van Allen belts 
is permanent as long as the satellite orbit passes, at 
least partially, inside them. 

We will show here that the charging is unacceptably 
high at otherwise favorable satellite orbits, so that some 
kind of charge removal technique is necessary, but this 
problem may be addressed rather easily by presently 
available technology. 

The range of the most favorable SEE orbital alti- 
tudes, roughly 1400 to 3300 km coincides with the 
inner region of the so-called inner radiation belt [21- 
24] , situated presumably near the plane of the magnetic 
equator. This region is characterized by a considerable 
flux of high-energy protons and electrons. For a SEE 
satellite at altitudes near 1500 km the duration of the 
charging periods is about 12 minutes. Maximum charg- 
ing rates occur in the central Atlantic. It should be 
noted that the South Atlantic Anomaly (SAA) — a re- 
gion of intense Van Allen activity which results from the 



low altitude of the Earth's magnetic field lines over the 
South Atlantic Ocean — cannot cause additional prob- 
lems for the SEE experiments. The reason is that the 
SAA mostly contains low-energy protons which cannot 
penetrate into the SEE capsule. 

Electrons are known to be stopped by even a thin 
metallic shell, so only protons are able to induce charges 
on the test bodies. Proton-induced charges on the test 
bodies can create considerable forces. The inner ra- 
diation belt contains protons with energies of 20 to 
800 MeV, and their maximum fluxes at an altitude 
of 3000 km over the equator are as great as about 
3T0®cm^^s^^ for energies E ^ 10^ eV and about 
2-104 cm-2s~i for E > 10^ eV. At 1500 km altitude 
these numbers are a few times smaller; the fluxes grad- 
ually decrease with growing latitude ip and actually 
vanish at ~ 40 ° . 

It is important that estimates of any resulting ef- 
fects take into account that (i) the capsule walls have a 
considerable thickness and stop the low-energy part of 
the proton flux and (ii) among the protons that pene- 
trate the capsule and hit the Particle, the most ener- 
getic ones, whose path in the Particle material is longer 
than the Particle diameter, fly it through and hit the 
capsule wall again. As for the Shepherd, its size is large 
enough to stop the overwhelming majority of protons 
which hit it. 

In what follows, we will assume a Shepherd radius 
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of 20 cm and a Particle radius of 2 cm and estimate the 
captured charges for some satelhte orbits in a capsule 
whose walls of aluminium are 2, 4, 6 and 8 cm thick. 
The SEE satellite must actually involve several coaxial 
cylinders for thermal-radiation control, and the com- 
bined thickness of their walls must amount to several 
cm. We will assume, in addition, that the Particle also 
consists of aluminium and stops all protons whose path 
is shorter than 4 cm (thus overestimating the charge 
by a small amount since most of protons will cover a 
smaller path through the Particle material). A 100 g 
Particle of aluminium will have a radius of « 2.07 cm. 

It is advisable to determine first which charges (and 
fluxes that create them) might be neglected. 

4.1. Admissible charges 

Let us estimate the Coulomb interaction both between 
the Shepherd and the Particle and between each test 
body and its image in the capsule walls. To estimate 
the spurious effects on the Particle trajectory, it is rea- 
sonable to calculate its possible displacements due to 
the Coulomb forces from the growing number of cap- 
tured charged particles. We assume that the test bod- 
ies are discharged by grounding to the capsule before 
launching the Particle in each given experiment. 

Criterion. We will call the induced charges, or the 
fields they create, admissible if they cause a displace- 
ment of the Particle with respect to the Shepherd 
smaller than a prescribed coordinate measurement er- 
ror 61 (we take here 61 = 10^® m) for a prescribed 
measurement time (we take t > 10* s). 

The Coulomb acceleration aqlt) ~ (lMq_rn/ {r'^fn) 
(in the Gaussian system of units) depends on the 
Shepherd-Particle separation r and on the form of the 
function J(t) , which in turn depends on the satellite 
orbital motion. 

The charge- induced Particle displacement is approx- 
imately 



a; = dt 



(19) 



since the acceleration is almost unidirectional. If, for 
estimation purposes, we suppose that the flux J(t, x) 
is time-independent, then the resulting displacement is 
about 



a; 



1 e^SAfSmJot'^ 

30 r^m 



(20) 



where Sm ~ 1256 cm^ is the Shepherd cross-section, 
Sm ~ 12.56 cm^ is the Particle cross-section, m is the 
Particle mass and r is an average Shepherd-Particle 
separation. 

The strong time dependence is explained by the 
rapid growth of the Coulomb force due to growing 
charges. Numerically, with the above values of Sm 



and Sm , taking m — 100 g and r = 1 m (the lat- 
ter leads to an overestimated force since the Particle 
spends most of time at greater distances), we find: 

JqI* < 0.83-1018 



2 —4 

s^cm ^. 



(21) 



For t ~ 10 s an admissible flux is thus less than 9 
cm s . 

Another undesired effect is that the Particle, being 
charged by the belt protons, will interact with the cap- 
sule walls. This is well approximated as an interaction 
with the Particle's mirror image in the wall, while the 
latter may be roughly imagined as a conducting plane. 
Then, assuming that the Particle is at average at about 
25 cm from the capsule wall and using the same kind 
of reasoning as above, we obtain instead of ( ^ ) 

J^t* < 2.07-1019 s^cm-^ (22) 

and an admissible proton flux less than 45 cm^^s^^ for 
t = 10* s. 

Some more estimates are of interest: if the charge 
can be kept smaller than a certain value, then what 
is the upper limit for it to create only negligible dis- 
placements? Suppose that there are constant charges 
on both the Shepherd {q = qu) and the Particle {q = 
q,n, rn — 100 g), then they are admissible according to 
the above criterion as long as 



qM^n 



< 2-10-6 CGSE/ = 1-10-24 
i < i-lO-^CGSE 2. 



C^ 



(23) 
(24) 



These inequalities follow, respectively, from considering 
the Shepherd-Particle interaction and the interaction 
between the Particle (located at 25 cm from the wall) 
and its image. Thus the maximum admissible Particle 
charge is about 7-10~* CGSE, sa 1.5-10^6; assuming 
this value, it follows from ( p^ ) that the maximum Shep- 
herd charge is about 3-10-3 CGSEg w 5.5-10^6. With 
these charge values the electric potentials on the test 
body surfaces are 



Um ~ 1.5-10-4 CGSE^/cm = 45 mV; 
[/,„ « 3.5-10-''CGSE,/cm = 105 mV. 



(25) 



If by any means the requirements (|2^) , ( |24| ) are sat- 
isfied (e.g., the potentials are kept smaller than the val- 



ues (25)), the electrostatic effect on the Particle trajec- 
tory may be neglected. 

The Shepherd's interaction with the respective im- 
age charge induced at nearest location to it in the SEE 
experimental chamber does not lead to appreciable Par- 
ticle displacements. A very demanding requirements on 
the Shepherd's charge emerges, however, if the SEE 
satellite is used for G-dot determination (whose de- 
tailed discussion is postponed to future papers). For 
that case, 



Um < 1 mV. 



(26) 
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Table 2. Average flux, peak flux and captured charges per revolution in some satellite orbits 



Orbit 


Wall 


Average flux, 


Peak flux. 


Shepherd 


Particle 






thickness 


— 2 —1 

cm s 


— 9 —1 

cm s 


charge qm 


charge q 


m 




2 cm 


1420 


12300 


e 


4.5-10^ 


e 


1500b 


4 cm 


1000 


8800 


MO^" e 


2.6-10^ 


e 




6 cm 


770 


6800 


7.5-10^ e 


1.5-10^ 


e 




8 cm 


600 


5400 


6-10^ e 


1.2-10^ 


e 




2 cm 


646 


5700 


6.5-10^ e 


1.9-10^ 


e 


1500c 


4 cm 


464 


4200 


4.3-10^ e 


l.MO^ 


e 




6 cm 


365 


3300 


3.5-109 e 


7-106 


e 




8 cm 


280 


2700 


2.7-109 e 


5-106 


e 



Evidently, in this case the Shepherd-Particle interaction 
per se is not the determining factor with respect to 
charge limits on the test bodies. 

4.2. The captured charges in certain orbits 

The charges captured by the Shepherd and the Parti- 
cle on board a satellite in various circular orbits for a 
single revolution around the Earth, a period of about 
two hours, were estimated in |2^. (Actual measure- 
ment times may exceed this period, but not by much.) 
These estimates were obtained with the aid of SEE2 
and SEREIS software, created at Nuclear Physics In- 
stitute of Moscow State University |2^ . 

The results obtained lead to some conclusions of im- 
portance for the SEE experiments (details see in p5| ). 

First, the models show zero proton fluxes in equato- 
rial orbits of 500 — 800 km altitude but indicate con- 
siderable fluxes at the same altitudes due to crossing 
the SAA. It turns out, however, that the SAA is over- 
whelmingly a low-energy phenomenon and leaves the 
fluxes virtually unaffected on the relevant energy scale, 
beginning at approximately 65 MeV. Moreover, there is 
a very small proton flux due to the SAA even at ener- 
gies above 10 MeV, hence, with 1 mm layer of shielding, 
the SAA influence is negligible. Behind a thicker layer 
of shielding there are actually no secondary particles 
due to SAA protons. 

Second, at an altitude of 1500 km the fluxes depend 
substantially on the orbit orientation but remain on the 
same scale of a few million protons per cm^ at energies 
over 65 MeV. 

Third, at an altitude of 3000 km both the total flux 
(for X = 0) and its high-energy component in particular 
are evidently a few times greater than at 1500 km. 

Fourth, and most important: for all orbits in the 
desirable range of altitudes the charges are quite large 
compared with their admissible values and they remain 
large even behind relatively thick walls. It is thus quite 
important to have means to detect and remove these 
charges during the measurements. Moreover, as seen 
from the peak values in Table 2 and time scans of Van 



Allen charging in orbits of interest (1500c is one of 
the most favourable ones, 1500b is less favourable; the 
data were also obtained using the above-mentioned soft- 
ware), at a charging peak when crossing the magnetic 
equator the time required for the charge on the test 
bodies to reach its maximum allowable values, as listed 
above, is a matter of seconds, not minutes. Therefore 
the charge must be detected and removed as it builds 
up, on a time scale of seconds. 

The detection and measurement of the charge on the 
test bodies can probably be achieved relatively easily 
by an array of minute microvoltmeters attached to the 
inner wall of the experimental chamber. 

Several methods for removing positive charge are 
now being evaluated. A simple and promising method 
may be to shoot electron beams directly at test bodies. 
The number of electrons needed is of the order of 10^/ s. 
Although this approach has the inherent drawback that 
it requires that an active system must perform correctly 
for many years, it is simple in principle and will accom- 
plish the goal. 

5. Concluding remarks 

The main results of the recent developments described 
in this paper may be summarized as follows: 

1. Numerical simulations of the Particle relative 
motion in the case of a lighter Shepherd (200 kg), for 
initial Particle-Shepherd separations between 2 and 5 
m, for Hoih = 500 3000 km, has shown that the sen- 
sitivity of trajectories with respect to changes of G did 
not change too much compared to the previous esti- 
mates with a heavier Shepherd (500 kg). In particular, 
variations 5G/G ~ 10"^ lead to trajectory shifts along 
the X axis ranging from 0.634-10"^ to 3.34-10~6 m. 

It has been found that an error of 1 cm in i?orb 
leads to trajectory shifts of 3-10^^ to 10^^ m, i.e., 
about an order of magnitude smaller than the planned 
coordinate measurement error, 10^^ jji^ and than the 
trajectory shifts due to G variations of 1 ppm. 

An estimated error 6G/G due to Shepherd non- 
sphericity, characterized by its quadrupole moment 
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J2 ~ 10^^, is close to 10^^; this implies rather se- 
vere requirements upon the Shepherd fabrication preci- 
sion. The same quadrupole moment causes even greater 
6G/G if the Particle motion starts at yo > 15 cm. 

2. Computer simulations have shown that the grav- 
itational constant G can be measured up to ^ 5-10^ . 
The inverse square law may be tested with a sensitivity 
of a ~ 2-10"^ for A = 1.2 m and a ~ 3-10"^ for 
A ^ 3.4-10^ m (half the Earth's radius). Observation 
of the Shepherd orbit precession makes it possible to 
test the inverse square law at A ~ 3.4-10^ m up to 
10-10. 

3. Estimation of test body charging due to crossing 
the van Allen radiation belts shows that this effect re- 
quires special means for charge measuring and removal. 
These means, however, do not go beyond the presently 

available technology level. 
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